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Motivated by the recent experiment [V.A. Sautenkov, Yu.V. Rostovtsev, and M.O. Scully, Phys. 
Rev. A 72, 065801 (2005)], we develop a theoretical model in which the field intensity fluctuations 
resulted from resonant interaction of a dense atomic medium with laser field having finite bandwidth. 
The intensity-intensity cross correlation between two circular polarized beams can be controlled 
by the applied external magnetic field. A smooth transition from perfect correlations to anti- 
correlations (at zero delay time) of the outgoing beams as a function of the magnetic field strength 
is observed. It provides us with the desired information about decoherence rate in, for example, 
^^Rb atomic vapor. 

PACS numbers: 32.80.Qk, 42.50.Ar 
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I. INTRODUCTION 

The fundamental limits of spectral resolution and 
sensitivity of spectroscopic techniques, the information 
transfer and computation rates, spatial resolution of op- 
tical microscopy and imaging are determined by sta- 
tistical properties of light. For the last five decades 
enormous theoretical and experimental research activi- 
ties have been devoted to studying fluctuations in classi- 
cal and quantum systems 1]. 

The first experiment on the correlation between the 
intensity fluctuations recorded at the two different photo- 
detectors illuminated by the same thermal light source 
was performed by Hanbury-Brown and Twiss ^ . In their 
experiment, photon bunching, i.e., an enhancement in the 
intensity-intensity correlations has been observed. 

Quantum formulation of optical coherences was intro- 
duced by Glauber in his pioneering work 0. Photon 
anti-bunching has been predicted by Carmichael et al. Q 
and then it was firstly observed in resonance fluorescence 
experiment by Kimble et al. j^. 

As a generalization of the results obtained for two- 
level atomic systems 0, to fluorescence from a A three- 
level atomic system showing an anti-bunching effect in 
the second order correlations has been studied in 8] . Due 
to four-wave mixing in cold atoms ^] under condition 
of electromagnetically induced transparency (EIT) [Toj . 
Harris and co-workers Tl^ have measured the correlation 
between Stokes and ant i- Stokes photons emitted from Rb 
atoms with short time delay. 

Kuzmich et al. 12] have demonstrated a generation of 
pair photons with controllable time delay in the issue of 
quantum information storage and retrieval 0, us- 
ing ensemble of atoms 15j. The correlated photons have 
been greatly attracted in the study of, e^, entanglement 
amplifier 16 1, subnatural spectroscopy |l3, quantum mi- 
croscopy 18j, nonclassical imaging of trapped ions 



and many others. 

A transition from anti-bunching to bunching of light 
emitted from a few atoms in a very high finesse cavity 
has been demonstrated [2^. The transition occur s by 
increasing a number of atoms interacting with light [2l|. 

The matched fields treated both classically and 
quantum mechanically 23] can be another promising 
theoretical approach to the switching of correlations in 
a three-level atomic sample. However, the result is 
very sensitive to the detuning between driven fields and 
atomic levels. Photon bunching in the intensity-intensity 
correlations between pump and probe fields for different 
probe detunings in Rb vapor (2^] and in the temporal cor- 
relations between forward and backward ant i- Stokes pho- 
tons scattered from sodium vapor has been demon- 
strated. 

The most recently, Scully and co-workers [2^ have ob- 
tained a smooth transition from EIT correlated to anti- 
correlated photons emitted from coherently prepared 
^^Rb vapor. 

In the present work we develop a theory to explain the 
results of the previous experiments [2^, whereas laser 
source to be considered here with a finite bandwidth. 
A diode laser used in our experiment, would have low 
intensity fluctuations but non-negligible phase fluctua- 
tions under certain condition. The fluctuations of the 
input light after interacting with the atomic sample can 
be enhanced and contain information about atomic sam- 
ple. For instance, this has been used as a spectroscopic 
tool [29|, |30|. Particularly, the laser phase fluctuations 
can be converted into the intensity fluctuations due to 
the interaction of field with atomic vapor 

Furthermore, based on the numerical results, we sug- 
gest a new promising method to estimate decoherence 
rate for Zeeman sub-levels. 

This paper is organized as follows. In the next section, 
the experimental setup is described and the obtianed re- 
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suits are reported. In section III, we study the absorption 
induced fluctuations of laser beam intensities and their 
correlations by considering a generic three-level A atomic 
system interacting with laser fields with orthogonal po- 
larizations. We obtain the approximate analytical solu- 
tions elucidating an origin of perfect correlations as well 
as anti-correlations between two modes. Then, in the 
next section, we solve numerically equations of motion 
to prove the results of the analytical predictions. The 
last section is conclusion. 



II. 



EXPERIMENTAL SETUP AND OBTAINED 
RESULTS 




PBS 



FIG. 1: A simplified schematics of experimental setup and a 
level scheme considered for Rb atoms. PBS: polarizing beam- 
splitter, A/2, A/4: wave plates, PD: photo-detector, DO: dig- 
ital oscilloscope, Comp: computer. 



A setup of the experiment (similar to one in 12^) is 
shown in Fig. 1. An external- cavity diode laser 6d\ is 
tuned to Di line (56'i/2(F = 2) ^ ^Pi/2{F' = 1)) of 
^^Rb. An input beam is separated by a beam-splitter. 
The polarizations of these two separated beams become 
orthogonal using a A/2-wave plate put on the way of one 
beam and these are combined together by a polarizing 
beam-splitter (PBS). After the A/4 wave-plate the beam 
is a combination of two circular polarized optical fields. 
A glass cell of length L = 7.5cm with Rb vapor (natu- 
ral abundance) at density approximately lO^^cm"^ is in- 
stalled in a two-layer magnetic shield. A simplified level 
scheme is depicted in inset of Fig. 1. The opposite circu- 
lar polarized beams interact with the vapor and induce 
a ground state Zeeman coherence in Rb atoms. EIT res- 
onance is presented in Fig. 2. Transmitted laser beams 
after the second A/4 wave-plate are separated again by 
another polarizing beam-splitter and focused on fast pho- 
todiods (PD) with frequency bandwidth 75kHz— 1.2GHz. 
The optical path lengths for both beams are the same. 
Signals from PDs are sent to a digital oscilloscope (DO). 
As varying a magnitude of longitudinal magnetic field 
the transmitted fields are changed at the optical power of 
0.5mW (total power ImW) and beam diameter of 0.1cm 
for each beam at the entrance window of the Rb cell. 



£ 0.50- 




B (Gauss) 

FIG. 2: A transmission of the optical field through the Rb 
cell as a function versus magnetic field B (solid curve). Dots 
stand for Lorentzian fit. The total optical power at all at 
entrance window is ImW. 



The time dependent intensity fiuctuations SIi^2{t) of both 




Time {\is) 

FIG. 3: Waveforms from photo-detectors with total optical 
power of laser beams at front window of Rb cell is ImW. 
The fluctuations of intensities versus time for two coinciding 
beams with (a) no magnetic fleld B = and (b) applied 
magnetic fleld B = — 0.47Gauss. 



optical beams transmitted through Rb vapor, are regis- 
tered by the photodetectors (see. Fig. 3). Data presented 
here is a part of the recorded data in 10/isec. The sig- 
nal in Volts is proportional to laser intensity as 500 V/W. 
Furthermore, the intensity-intensity correlations between 
two modes can be calculated using the observed data for 
the intensity fiuctuations. The second order correlation 
function G^^^ (r) for intensity fiuctuations of two optical 
beams with time delay r is given by 



G(2)(r) = 



{Sh{t)Sl2{t + T)) 



(1) 



where the time average of arbitrary variable Q{t) is de- 



fined as {Q{t)) — J^~^^ Q(t)dt/T. The integration time T 
is taken to be as large as 10/is. In the absence of the exter- 
nal magnetic field 5 = where EIT condition is fulfilled 
(two-photon detuning is zero), the induced fluctuations 
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of the transmitted beams by Rb vapor are almost syn- 
chronized (see Fig. 3(a)). In the case of zero detuning, 
the intensity- intensity correlation curve of Fig. 4(a) has 
a sharp spike clearly showing bunching. The magnitude 
of the correlation peak at r = is of 0.9 and the aver- 
age background is near 0.15. The width of the correlation 
peak increased as reduces the optical power. On the other 
hand, the most intriguing feature is observed when an 
applied magnetic field is of 5 = — 0.47Gauss. As is seen 
from Fig. 3(b), intensity fluctuations are out of phase 
when two photon detuning becomes non zero. A presence 
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FIG. 4: Intensity-intensity correlation functions G'''^\r) as 
functions of time delay r for magnetic fields B = (above) 
and B = -0.47 (bottom). 



of the magnetic field demonstrates an exhibition of anti- 
correlation (see Fig. 4(b)) as is expected from the data 
plotted in Fig. 3(b). The width of the peaks is associated 
with the saturated width of resonance in Rb vapor ab- 
sorption (a single photon resonance) 36]. Moreover, the 
second order correlation G*^^^(r) obtained with spatially 
separated beams (distance between beams 0.3cm which is 
bigger than the beam diameter 0.1cm) has a correlation 
peak (at r = 0) of 0.7 and a larger background of 0.3. 
This contrast indicates that the enhanced correlations are 
clearly due to overlap of two beams. We have also per- 
formed a set of measurements of G^^^ (r = 0) for different 
values of the magnetic field B at optical power ImW (see. 
Fig. 10). A measurement shows that a perfect switch- 
ing from photon correlation to anticorrelation where the 
correlation peak with magnitude 0.9 for zero magnetic 
field B = switches to a dip with magnitude —0.9 for 
B = — 0.47Gauss. This pronounced modification of wave- 
forms is very important in determination of some experi- 
mental parameters which will be shown later. The width 
related to intensity- intensity correlations G*^^^ (r = 0) are 
0.24Gauss, almost four times narrower than the corre- 
sponding EIT width which is 0.85Gauss for optical power 
ImW. 



III. THEORETICAL RESULTS 

A. Three-level A atoms driven by laser with finite 
bandwidth 

Let us consider interaction of two modes generated by 
the diode laser with three-level A atoms that have a level 
scheme shown in Fig. 1. The equations of motion for this 
system are given by (e.g. see 8j) 

Pbc = -{-f3^i^)pbc^i^lpla-^^2pba 
Pba = -(71 + ^Ai)p5a - i^liPbb - Paa) 

-iVt2pbc + iHt)Pba 
Pea = -{l2-i^l)Pca-i^2{Pcc- Paa) 

-iQipl^^i(j){t)pca 

Pbb = IsPcc^llPaa^i^lpla-'^^lPba (2) 
Pec = l2paa - IsPcc - i^2pca + ^^2pL 

where paa = ^ — Pbb — Pec The single and two photon de- 
tunings are Ai = uJa—uJc—J^2 = i^i— c^^+c^b, (z^i = 1^2) and 
A = uJb — LJc' The effective decay parameters defined 
as 73 = (73 + 721 + 7i2)/2, 71 = (71 + 72 + 721 + 73i)/2 
and 72 = (71 +72+73 + 712+732)72, where 71 and 72 
correspond to the spontaneous emission rates from level 
I a) to levels |c) and respectively; 721, 712, 732 and 
731 stand for dephasing rates, and 73 is population de- 
cay rate of the level \b). In derivations of Eq.(|2|), we 
have kept the operator normal ordering, i.e., we neglect 
rapidly oscillating terms 6] and used re-scaled variables 
as pij 1/Ncpij^ Nc is number of collective atoms. A 
diode laser radiation experiences phase diffusion, and the 
phase <j){t) in Eq.Q represents the fluctuating phase of 
driven field which is characterized by Wiener-Levy diffu- 
sion process ^] . For such process average and two-time 
correlation function of stochastic variables are given by 

iW) ) = 

m)<P{t')) = 2DS{t-t') (3) 

where D is the diffusion coefficient; the stochastic aver- 
ages denoted by the upper bar. Thus, the input laser field 
has a Lorentzian spectrum with a FWHM bandwidth of 
D/ttHz. In a realistic situation, the phase correlation 
has a finite relaxation time. The Gaussian process in 
which the correlations are determined by the exponential 
function of time delay is often referred to as Ornstein- 
Uhlenbeck ^] or colored noise. A stationary equation 
after taking stochastic average of Eq.(|2|) is shown in Ap- 
pendix. The numerical simulations of Eq.(|2|) will be pre- 
sented below. 
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B. Absorption induced intensity-intensity 
correlations 



Propagation equations for the laser fields are given by 



dz 



dz 



(4) 



In order to give a qualitative theoretical analysis of 
our experimental results, let us adopt a theory which im- 
plies for a thin absorbing medium. It is assumed that the 
transmitted field could be understood as a superposition 
of input and induced fields in the first order approxima- 
tion for /^i,2^, if 1^1,2^ <C 1; here L is the length of the 
atomic sample and i<ii^2 are some coefficients 30j. Fur- 
thermore, this could be still valid for a preferably long 
medium with a weak absorption, but satisfying the con- 
dition /^i,2^ <C 1. In what follows, we will show that 
this approximation reproduces the observed results qual- 
itatively. Following Walser et al. 30] and Martinelli et 
al. the transmitted fields are given by 



Qi -^i{niL)pac{t), 

9.2 ^i{K2L)pab{t)' 



(5) 



To include Doppler effect, the coherence terms in Eq.(|5|) 
should be averaged by the Maxwell-Boltzmann veloc- 
ity distribution. However, the Doppler broadening may 
play important role in many other experiment with Rb 
atomic vapor, but, in this case, it turns out not to be 
so crucial, because we are interested only in correla- 
tion behaviour as functions of two-photon detuning, in- 
stead of one-photon detuning. The transmitted intensi- 
ties Ii,2{t) ^ |^i^2^(^)P are, thus, written as 



hit) 
hit) 



VLiiKiL)lm{pacit)} 
Q.2in2L)lm{pahit)] . 



(6) 



here paq is the atomic coherence term for level a and 
{q = 6, c) and we assume that the input fields are real 
and much stro nger than the induced ones. Defining that 
6Q{t) = Q{t) — Q{t) stands for the fluctuation of arbitrary 
variable the intensity fluctuation to be read 



5hit) = ni{KiL)lm{Spacm 
Shit) = Q2{^2L)lm{5pabm. 



(7) 



From this expression, it is easy to check that {SIi^2{t)) = 
0, because {Spaq{t)) = 0, g = c, 6. From Eq.Q, using 
the definition Eq.Q one obtains the intensity-intensity 
correlations 



G(2)(r) 



{lm{Spaeit)}lm{Spabit + t)}) 

^{[lm{6pacit)W){[lm{Spabit)W) 



(8) 



The undertaking process is stationary, thus, the argu- 
ment t+r of the last term in the denominator is displaced 
by t. 



C. Approximate theoretical analysis 

Zero detuning 

Let us consider the resonant case where all detunings 
are set to be zero and l^i = (^2, 73 = 0, 71 =72- It is easy 
to show analytically that (p^^) = (p^^) where an equation 
for stationary state is given in Appendix. According 
to exact numerical simulations of Eq.(|2|), the coherence 
term p^c are real i.e., p^c = pj^- Therefore, equations for 
coherence terms p^a and pea can have symmetrical forms 
as 



Pba = -llpha - i9i{pbb - paa) + i4>{t)Pba - ^^iPbc 
Pea = -1lPca-i^l{pcc- Paa)^iHt)pca-i^lPbc 

(9) 

From Eq.Q, one can see that two equations, thus, two 
modes are decoupled. Note that because of symmetrical 
atomic configuration, it is obvious that populations p^b 
and pec are identical P55 ~ pec This is, of course, true 
only in a resonant case. From Eq.Q, it follows that 



Pba{t) ^ Pca{t) 



(10) 



This correlated behavior could be understood as follows. 
Phase fluctuations of incident beams are converted into 
intensities fluctuations via atom-field interactions as is 
seen from Eq.Q. Roughly speaking, the three- level 
atoms would experience driven fields with the effective 
Rabi frequencies fluctuating around which is resonant 
to the degenerate lower levels. Because the noise contri- 
bution is the same in two modes, any instant deviations 
from the resonance condition will be also the same. Thus, 
the induced absorption should be also the same for both 
modes. In this sense, this system would be in close rela- 
tion to what is called a correlated emission laser firstly 
proposed by Scully 1], in which pairs of induced pho- 
tons of different modes can be generated simultaneously 
exhibiting a sharp bunching. 



Non-zero detuning 

In the non-degenerate situation, the equations for two 
modes are coupled. However, the effective Rabi frequen- 
cies would have again the same fluctuations, but, Qi is no 
longer resonant to the lower levels. The induced absorp- 
tion should not be the equivalent in this case, because, 
the deviations of the effective Rabi frequencies would be- 
come farther from one of ground levels, but, closer to 
another at any instant time. As a consequence, the pop- 
ulations of the excited and ground states would also fluc- 
tuate. By virtue of P55 + pee ^ 0, since paa ^ 0, it is 
possible to do the following assumption as 



Pbb{t) 

Pcc{t) 



■fit), 
fit)- 



(11) 
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where a fluctuation f{t) is real and c is constant. Note 
that it is not necessary to know an exphcit expression for 
f{t). This couphng function f{t) appears only because 
of non-zero detuning, otherwise it is zero. The equations 
for two polarizations, p^a and pca, are given by 

Pba = -{jl^iAi)pba-i^lf{t) ^i^{t)pba-i^lPbc^c' 
pea = -{jl-iAi)pca^i^lf{t)^i^{t)pca-i^lplc^c' 

(12) 

The formal solutions to be read 

Pba{t) - +<^>(t)+Ci 

Pca{t) ^(t)+C2 (13) 

where ^{t) = -zl^i /^'^ dt'e-^i(^-^')+^(^(^)-^(^'))/(tO, 

Ci,2 = -^^^1 Jl 6^t'e-^i(^-^')+^(^(^)-^(^'))(p5c,c6 + C) which 

can be slowly varying for to ^ — oo, t ^ oo. Note that, 
contributions of Ai to solutions are neglected since, later 
on, only imaginary part of amplitudes will be of interest. 
Ea. (|T3j) clearly indicates an exhibition of anti-correlation 
between two modes. 



D. Numerical results 

In the Ornstein-Uhlenbeck process j^, the colored 
noise (_{t) yields the steady-state correlation function 

(mm) = eAie-^-l*-*'l (14) 

with (^(t)) = 0. The stochastic differential equa- 
tion Eq.© can be solved using Monte-Carlo numerical 
simulations. A Box-Mueller algorithm and the Euler- 
Maruyama method are used to realize the colored noise. 
Namely, one can see that the generated noise by the fast, 
integral algorithm developed in , is in a perfect agree- 
ment with the analytical deflnition given by Eq. (|T^ due 
to averaging over as many as 1000 different realizations. 
A relaxation time 1/Xl is taken to be small, to ensure 
that undertaking process would be approximately white 
noise, i.e., ^(t) ~ ^{t). As a matter of fact, various 
choices of parameters © and Al, should not drastically 
influence to the flnal results. The numerical solutions 
of Eq.(|2|) allows one to obtain the intensity fluctuations 
deflned by Eq.©. The numerical results of Eq.© are 
plotted by dotted and solid curves in Fig. 5. In resonant 
case, dynamics of two modes are in phase i.e., well syn- 
chronized. The dephasing rate for both cases are chosen 
to be much smaller than decay rates. Absolute values of 
two Rabi frequencies are the same as is considered in the 
experiment. If two-photon detuning A becomes compa- 
rable to 73 then dynamical behaviors are absolutely out 
of phase. It is seen more clearly from Fig. 6. Dynami- 
cal features shown in Fig. 5, can be seen more clearly in 
terms of cross correlation functions. Figure 6 describes 
a switching between two completely different behaviors, 
namely, correlation and anti-correlation of two modes. 
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FIG. 5: Synchronized dynamics for intensity differences of two 
fields (dotted and solid curves) given by Eq.(|2| are depicted 
for resonant (above figure, A = 0) and non-resonant (bottom 
figure, A = O.OI71) cases. All rates and Rabi frequencies are 
taken to be 712,21 = O.OI71, 713,31 = 73 = 0, 71 = 72 and 
Qi = Q2 = 71 
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FIG. 6: A perfect photon correlation (a spike at r = 0, see, 
(a)) for A = and an ant i- correlation of two modes (a dip 
at r = 0, see, (b)) between two modes for A = O.OI71 are 
obtained. The plots are calculated for the same parameters 
as in Fig. 5. 



IV. A DECOHERENCE RATE DETERMINED 
BY SWITCHING IN INTENSITY-INTENSITY 
CORRELATIONS 

In what follows, we analyze this switching in more de- 
tail. Similarly as in 28], we focus on correlation functions 
with zero time delay G^'^\r = 0). First of ah, Eq.Q 
is obtained for fixed fti but, different 73 and shown in 
Fig. 7. Note that transitions from correlated photons to 
anti-correlated ones are appeared to be smooth and have 
certain widths. These widths are getting more wide with 
the increase of dephasing rates. Surprisingly, as shown in 
Fig. 8, the widths are 'invariant' as functions of re-scaled 
detuning variable A/73 with respect to the correspond- 
ing decoherence rates 73. Moreover, let us test also how a 
Rabi frequency's change might affect to correlations. In 
Fig. 9, we depict numerical results for Eq.(|S|) depending 
on not only detuning, but also, Rabi frequencies for two 
different fixed values of 73. As a matter of fact, the corre- 



-0.6 -0.4 -0.2 0.2 0.4 0.6 

B (Gauss) 



FIG. 7: 3D plot for correlation G^'^\t = 0, A,73) with zero 
time delay (r = 0) as a function of detuning A and decoher- 
ence rate 73. Rabi frequencies are Qi = Q2 = 71. 
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FIG. 8: A sample of correlation functions G^'^\t — 0, A/73) 
with scaled variable as A/73 for different decoherence rates: 
73 = 0.01, 0.02, O.II71. All curves coincide. The plots are 
calculated for the parameters taken from Fig. 6. 
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FIG. 9: 2D view of 3D plot for correlation G^'^\t = 0, A, Qi) 
depending on different values of both detuning and Rabi fre- 
quency for fixed dephasing rate 73 = O.I71. 



lation curves are again 'invariant' for all Rabi frequencies 
those being not smaller Vti > 71. These two intriguing re- 
sults for the correlation functional invariance do promise 
a relatively precise determination of the decoherence rate 
in the robust way from the experimental point of view. 
As is mentioned in the introduction, the atomic energy 



FIG. 10: A finding dephasing rate in Rb vapor. Correla- 
tion G^'^\r = 0,B) varying with magnetic field strength B 
is compared to the experimental data which allows us to get 
scaling factor being a ~ 1[G']. Dephasing rate is, thus, found 
to be 73 ^ IMHz independently from 71. No free parameters 
are used here. Experimental data is presented by circles for 
optical power ImW 



levels are perturbed, due to the interaction of the mag- 
netic moments of atoms with the external magnetic field 
B. This leads to a non-degeneracy of atomic ground lev- 
els. This shift is given by 

A = aB (15) 

The constant is defined by Bohr's magneton, the mag- 
netic quantum numbers and gyromagnetic factor as a = 
{iJiB /K)g{m2 — mi). To be more explicit, let us con- 
centrate on ^^Rb vapor in connection to our experi- 
ment [2^. An external- cavity diode laser is tuned to 
^Si/2{F = 2) ^ 5Pi/2(F' = 1). Applied magnetic field 
leads approximately to three level A atomic configura- 
tion with a common upper level 5Pi/2{F' = 1,^3 = 
0) and two lower levels ^Si/2{F = 2, mi = —1) and 
bSi/2{F = 2,m2 = +1). Using the facts that g = 0.5 for 
Si/2{F = 2) and iib/^= 1.4MHz, the constant would be 
estimated a = 1.4MHz/G. Moreover, from Eq.(^3, the 
relation for B can be rewritten as 

A 

B = a— (16) 

73 

where a is a scaling factor. For fixed a, all correlations 
with no time delay are supposed to be invariant. As is 
seen from Fig. 10, experimental data (circles) are form 
preserved and have the identical widths regardless of op- 
tical powers choices. Asymmetry with respect to the zero 
detuning may be due to Stark shift which is not of in- 
terest in the present situation. Because of two condi- 
tions of invariance, there is a very good reason to believe 
that a form of theoretical correlation functions for all 
^1 ^ 71 and any of 73, should be equivalent to the ex- 
perimental data measured for particular ftexp and jsexp- 
In this spirit, theoretical curves by changing scale a, can 
be compared with the measurements. Remember that 
= Isexp/a^ a new formula for dephasing rate can be 
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given as 

73exp = C^Ci (17) 

Once a could be found from experimental data, the de- 
phasing rate for that atomic vapor, is determined by for- 
mula Eq. ([T7|) . As we expected this is also independent 
from population decay rates. For Rb vapor, we have 
found from Fig. 10, that scaling factor a ^ 1[G]. There- 
fore, the decoherence rate is estimated to be jsexp ^ 
1 MHz. 



V. CONCLUSION 

An experimental demonstration of intensity correla- 
tions and anti-correlations of coupled fields in a dense 
Rb vapor is reported. A lower level coherence is cre- 
ated between Zeeman sub-levels by two laser beams with 
orthogonal circular polarizations. Intensity fiuct nations 
induced by resonant medium are correlated under reso- 
nance EIT condition and anti-correlated in presence of 
non-zero two photon detuning. A narrow correlation 
peak and anti-correlation dip, in time domain, are as- 
sociated with frequencies above EIT width and natural 
optical width. A dependence of correlations on mag- 
netic field (two-photon detuning) show resonance behav- 
ior. The resonances are near 4 times narrower than the 
width of the observed EIT resonances. A smooth transi- 
tion from perfect correlations to anti-correlations (at zero 
delay time) between the outgoing beams as functions of 
the magnetic field strength is robust with respect to a va- 
riety of different choices of physical parameters involved 
and, thus, can provide us with the desired information 
about decoherence in three level atomic vapor. More- 
over, correlation properties of coupled fields in A scheme 
can be used to reduce noise and improve performance of 
EIT based atomic clocks and magnetometers. The phase 
noise to intensity noise conversion is an important phys- 
ical process limiting the accuracy. In EIT atomic clock 
and magnetometers j42l 0, 10, it is possible to avoid 
the contribution of the atomic medium induced excess 
intensity noise. 
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Appendix: Stochastic averaging of equations with 
multiplicative noise 

A stochastic average of an arbitrary dynamical variable 
F(x), i.e., a path-integral over all possible realizations of 
the random numbers x(t), is given by 

A white noise does satisfy the relations x{t) = and 
x{t)x{s) = TS{t — s). A functional measure Djii[x] has a 
Gaussian density 

Di^[x] = N-^Dxe-^^ 
here N is normalization coefficient to assure 

I D^i[x] = 1 

A characteristic function is found to be 

Z{g) = e^IdrxiT)giT) = ^ (18) 

Following Wodkiewicz let us consider the following 
stochastic equations 

— Mo^ + ix(t)M^f 

In the interaction picture where ^ i{t) = e~^°*^(t), the 
equation can be rewritten as 

= ix{t)Mj{t)^j (19) 

where M/(t) = e"^o^Me^o^ A formal solution of 
Eq. (|T^ is given by 

here T is the time ordering operator. Using the relation 
Eq. ([T^ . the stochastic average to be read 

This is equivalent to the equation 

^5i = -fM,W'*, (20) 

and, finally, we arrive at 

^ = Mo*-|m2*. (21) 

This is the expected stationary equation. 
* e-mail: goa@physics.tamu.edu 
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